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Summary 


The flow over wave shaped surfaces is studied using a Navier Stokes 
solver. Detailed comparisons with theoretical results are presented, 
including the stability of a laminar flow over wavy surfaces. Drag charac- 
teristics of nonplanar surfaces is predicted using the Navier-Stokes solver. 
Also discussed are the secondary instabilities of wall bounded and free 
shear flows. 
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Section I 


Introduction 

At NASA/LaRC there has been an ongoing research program related to flow 
over wave shaped geometries. The primary objective of this program is to 
develop surface geometries that have lower aerodynamic drag compared to 
planar surfaces. This report is an outgrowth of an ongoing computer 
simulation study of flow over wavy surfaces. The details of the theoretical 
model and test results of the simulations for both laminar and turbulent 
flows have already been documented in two previous contractor reports (Ref. 

1, 2; NASA CR159305, 1980; NASA CR3408, 1981). Extensive simulation studies 
that were conducted since then are the topic of the present study. 

The problem of flow over wavy surfaces has been a major area of 
research in fluid mechanics for over a century. The mechanism of wave generation 
and growth of small amplitude waves in the sea (air-sea interact ion) has been a 
subject of sustained study for many decades. The appearance of ridges in sand dunes 
is of interest in sediment transport technology ( porous wavy surfaces ) . Small 

amplitude waviness of machined surfaces and their effect on the performance 
of aircraft (LFC aircraft for instance) is another application of 
wavy wall flow phenomena ( transition ) . For sinusoidal wavy surfaces in a 
narrow band (a/A ~ 0.03, 6 /A ~ 0(1)) it has been observed experimentally 
that the wave average skin friction drag in turbulent flow is substantially 
lower than that of a flat surface. However, the associated wave drag penalty 
for sine surfaces more than over-compensates for this drag reduction, 
resulting in a net drag increase for these surfaces. With the availablitiy 
of high speed computers it has been possible in recent years to solve many 
intractable problems. Rather than a hit-miss experimental study for low drag 



wavy surfaces a reliable computer simulation program (which is economical 
and accurate) was developed to analyse the total drag properties of 
both symmetric and asymmetric waveforms. Details- of the computer 
algorithm (Navier-Stokes solver) are given in Ref. 1. Application to 
turbulent flow and comparison with wind tunnel data are given in Ref. 2. 

The present work examines in detail the efficacy of the Miles-Ben jamin 
theoretical model (3) for laminar flows. In Sec. 2, comparisons of laminar 
flow calculations with theoretical models are presented. In the presence of 
disturbances, a steady laminar flow undergoes transition at a certain 
critical Reynolds number. The presence of waviness on the surface can itself 
produce these disturbances when the wavelength of the surfaces (a ) is in the 
critical range (Sec. 3 ). For the computer simulation to be effective in 
rating various geometries (based on their total drag) a proper comparision 
of methodology is required. Section 4 discusses this concept in detail. 

In Sec. 5, the results of laminar flow' over some low drag surfaces are 
discussed. The results of drag computations for turbulent flow over a 
variety of asymmetric and symmetric surfaces are discussed briefly in Sec. 

6 with the secondary instability of wall bounded shear flows considered in 
Sec. 7 . A brief discussion of the work carried out in the development 
of compressible SALLY program under the present contract and work on 
free shear flows is reported in Sec. 8. 
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Section II 


Laminar Flow Over Wavy Surfaces 
Theoretical Models 

The simplest model of flow past a wave train is based on the classical 
Kel vi n-Helmholtz theory (4). In this analysis, viscosity is neglected, 

and the primary flow Is uniform. Therefore the disturbance (waviness) 
can be studied through potential theory. The perturbation pressure on the 
wave is proportional to the square of the relative velocity between the wave 
train and the fluid flowing over it and is in antiphase with the wave 
elevation. Even though this model Is naive physically, the pressure component 
predicted by Ke I vi n-He Imholtz theory Is consistent with the results of a more 
realistic theory (3) valid for thin boundary layers, (Re -+°° or 6/X<<l.) 
Benjamin (3) shows that this pressure component is the only stress remaining 
in the limit as Re + « uniformly. 

Miles ( 5 ) presented an analysis of the wind wave problem considering the 
inviscid perturbati ons to a primary profile due to the disturbances. The 
primary profile has the effects of viscosity (profile curvature) in this 
analysis. For travelling wave forms where the critical point is away from 
the boundary a substantial component of pressure in phase with the wave slope 
(sheltering) may occur. Kendall (6) has experimentally obtained qualitative 
agreement with Miles theory. 

Benjamin's Theoretical Analysis 

Benjamin (3) developed a uniformly valid asymptotic analysis of the flow 
past a wavy boundary. Benjamin's analysis is based on the following 
assumptions (i ) the amplitude 'a' of the wave is much smaller than the wave 
number k (=2tt/a, where x = wavelength), (ii) the primary flow is parallel; 


3 


i.e. the Reynolds number Re^ is large so that the boundary layer 
growth over several wavelengths is miniscule, i.e., U(r|,£) = U(r|) 

where ri,£ are the normal and tangential coordinates of a point. 

Ihe total velocity field is assumed to be derivable from a stream 
function, i.e., 

¥(€» n) = ^ o (n ) + a{F( n ) + (U(n ) - c) e~^} e ik £ 
where 4 , 0 (n ) = jo 1 ^ LU(n ) - cj dp and c is the wave speed of the surface. 

The perturbation due to the wave, F(n ) is influenced by viscosity in a 
friction layer whose thickness can be evaluated for given flow conditions. 
This analysis yields expressions for perturbation stresses (pressure 

and shear stresses) for both laminar and turbulent flow over wavy 

, 3 

surfaces. 

Since Benjamin’s work a large number of experiments have been 
conducted over wavy surfaces (turbulent flows). Discrepancies noted 
between the experimental measurements and theoretical results prompted 
several workers to develop theoretical models which include perturbation 
Reynolds stresses (7,8) and solution of the full Nav i er-Stokes equations 
(9, 10, II, 12, 13). 

There are very few (reported) laminar flow experiments over 
wavy surfaces. However, laminar flow calculations serve to calibrate 
both theoretical and computational methods perhaps better than the 
turbulent flow case chiefly because of the turbulence closure problem. 
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Therefore we have performed a number of laminar calculations and 
compared them to the threoretical model of Benjamin (3). 
Navier-Stokes Solver 


Our Navier-Stokes solver utilizes conformal mapping and periodic 
boundary conditions in the flow direction. An inverse stretching 

extends the computational lid in the normal direction. Under these 
conditions the unsteady solutions obtained using this code 

represent the proper physical evolution of the flow. A decided advantage 
of the periodic code is that, compared to an inflow-outflow code, the 
periodic code is more economic for problems related to wave trains. 

Using a Rayleigh analogy concept, which is discussed in Sec. 4, the 
periodic code has been utilized to study total drag character! st ics 
of various wavy surfaces. Comparison with experiments are presented 
in Ref. 2. The agreement between the Navier-Stokes simulations and 
wind tunnel experiments has been satisfactory for the range of 
geometric parameters and flow parameters tested thus far. 

The two dimensional time dependent Navier-Stokes equations are: 

p (5 u/3 t + u • V u ) = -V p + V • x (2.1) 

V . u = 0 (2.2) 
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where t-jj = cr-jj, (which includes both laminar and turbulent stresses for 
turbulent flows). 

The domain of analysis and the appropriate boundary conditions are D : 

0 < x < 2r , f(x,t ) _< y ^ 00 

where the wavy surface is 

y w = f (x.t) 

Periodic boundary conditions apply in the flow direction. The boundary 
conditions appropriate for the normal direction are 

U(x ,o ,t ) = Uw(x,t) 
v(x,o,t) = V w (x , t ) 

and U(x,«,t) = V(x,«,t) = p(x,°°,t) = 0 
for a stationary wall 

Uw = V w = 0 

For a wall with normal suction velocity V w and are nonzero, related to 
the components of the suction velocity i.n the physical x, y coordinate 
direcions. Similarly nonzero components of velocity at the wall arise when 
the wall is moving with a phase speed. 

The domain bounded by the wavy surface (5) is mapped to a rectangular 
E -n coordinate system using conformal mapping techniques (1) i.e., to the 
domain Dj : 

0<X<3t, -1<Yj<+1 

Initial Conditions 

At time t = t 0 we prescribe a velocity field appropriate to a starting 
Reynolds number at the beginning of the wave section., i.e., 
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U(X,Y,t 0 ) = U 0 (Y) (2.3) 

The solution algorithm then develops a quasisteady solution appropriate 
to the flow over wavy surface. 

for laminar flows U 0 (Y) will correspond to the Blasius profile 
appropriate for the starting length. For turbulent flows, a starting profile 
related to the virtual length of the flow is prescribed. In later sections, 

(3, 4) the properties of the time evolving flow solutions are related to 
relevant physical phenomena. For the present, we briefly state that whenever 
a stable solution for wavy wall flow exists, the solution obtained under 
these initial conditions and boundary conditions evolved In a fashion similar 
to the solution over a wave traln(wlth a mean growth of boundary layer thickness 
due to the downstream length). For Instance, the pressure distribution over the 
first wave will not be quite the same as that a few wavelengths downstream because 
the flow has a mean boundary layer growth and therefore the surface has a 
smaller (x/S) or (a/5). The quasisteady solutions of the Navler-Stokes 
equation show these variations with time. In section 4, we present a detailed 
discussion of the evolution In time of the periodic code described herein 
and evolution of flow over length for an Inflow-outflow code over a 
wavetraln. 

For the test problems that are reported herein the results for pressure 
and skin friction show minimal growth with time mainly because the 
wavelengths considered are quite small compared to the starting length (x/X ~ 
0(10)). Also, the time of evolution where quasi steady state is reached 
(related to x in an inflow-outflow code) is quite small for these cases and 
hence comparisons with Benjamin’s theory are appropriate. 
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Critical Examination of Benjamin's Theoretical Result 

Benjamin develops results for rigid walls (c = 0) based on three 
special cases: 

(a) Assumption of a I i near profile U(Y) = GY; 


The perturbation stresses are given in (5.6) and (5.7) in Ref. (3) and 
are written down in nondimensi onal form (after algebraic manipulations) for a 
laminar boundary layer as, 
pressure : 

C p = {2P s /p U 2 ^ ) = -.021308 kaRl/2( X/> )-2/3 fcos(kx .„ /6 ) (2.4) 


+ P cos (kx - 7 t / 3 )} 


shear stress: 


Cf = (2t S /p l^oo ) = 


0.143806 ka( X/A )“^/3 { cos (k x +tt/6) +P cos kx) 
p. = 3.818 (X/A ) 2/3 R-i/2 


(2.5) 


Admittedly, the linear profile suffers from a lack of realism. 


However, when the wavelength of the wall Is extremely short, the disturbance 
would be expected to penetrate only a very small distance away from the wall, 
where a linear approximation to the velocity profile is justified. 

(b) A sinusoidal profile approximating the Blasius profile; i.e., 

U = sin KY; 0 < Y _< 5 (2.6) 

U = ; Y > 5 (2.7) 

K = (tt / 26 ) , 6 = 5 X R" 1 / 2 (2.8) 


Benjamin derives expressions for pressure and shear stress, for the case 
K > k i.e., 5/A < 1/4; the expressions are. 
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Cp = R % [(2aK 2 /* X (1 + k /i) x e 1 ( kx + " )] 

L(1 - k/£ cot 0 ) + 1.288 e^/fyn j% x (1 + k/£ ) ] 

C f = Rfc [ 2a K 2 k/£ m ei( kx + * 3 ) x (1 + k/£ )] 

[(1 - k/£ cot e ) + 1.288 x (1 + k/£)] 

where m = (kKUo/v)l/ 3 and 6 = £6 = ir/2 /I - (k/K) 2 

When 6 /x > 1/4 the analysis presented by Benjamin can be suitably modified to 
obtain the expressions for and as 

C p = Rfc [2K 2 a/jj, ( ( 1+n/l-n ) + e~3o ) e 1 "^ + O] 

[{ (1+n/l-n) (1.288e i7T / 6 m/£ - 1) + e" 25 (1 + 1. 288e ilT / 6 m/t )} ] 

C f = R£ 2.746 i^ak/jim ((1+n/l-n) + e-^ ) e 1 ( kx + * / 3 ) 

[{ (1+n/l-n) (1.2886^ /®m/£ - 1) + e“^ (1 + 1.2886^ /®m/£ )} ] 

where 0 = £6 = tt/2 / (k/K) 2 - 1 ; n = k/£ and m = (kKl^/v )1/ 3 ** 

(c) For the case of a Blasius profile or any general profile, when 
A = (kml/oo/U* (o) 2 ) = (km/K 2 ) € 1, Benjamin develops expressions for Cp 
and Cf as 

Cp = Re, L2k 2 a ei( kx +0(1 - 1.28& e 1 */ 6 ) U 2 e - kY dYJ 
Cf = Ffc [2.744 k 2 a/m ei( kx + tt /6) j 

**For a Blasius profile, 

K = 0.332 X Rj!/ 2, m = 0.20336kRl/ 2 (X/x )’ 2/3 

£ = 0.332 Rl/2/x ✓ (358(X/X ^R" 1 ) - 1, 

and 0 = 1.57 / 358(X/X ) 2 R" 1 - 1 
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The Implication of these various approximations Is now 
examl ned. 

(i) Assumption of small "ka" limits the accuracy of the solution to small 
amplitudes. 

Benjamin suggests ka should be - 0.02. 

(ii) Assumption of parallel flow requires to be large I.e. — to be of 

Y 

the order of 10. (For the 9 rowth of boundary layer over a wavelength is 

of the order of 5%. ) 

CHI) Of the three sets of expression that Benjamin develops for a Blasius 
flow, case (b) seems to be the most appropriate one. Case (a) holds 
only for very short wavelengths. Case (c) lacks realism as will be 
shown below. Case (b), where a sinusoidal profile is used to approximate the 
Blasius profile is a more realistic approximation in that the differences between 
a sinusoidal profile and a Blasius profile with the same wall slope 
(k = 7 r/ 26 ; 6=4.73 /vx/U^) are only minimal and should not cause an appreciable 
difference in the outcome of the results. 

The results obtained using case (c) may themselves be of little use. Here 
the basic requirement for the results to be valid is that the parameters 

m6 = 4.2 (KX Q ) 1/3 » 1 (2.17) 

and A = 6.3 (KXq) 4 ^ 3 R' 1/2 « 1 (2.18) 

The first of these conditions must per- force be met since assumption of oarallel 
flow requires X Q /A » 1 and hence kX Q = 2ir ( k = 2 tt/A ) . However the requirement 
that A « 1 may not be met for all flow condition. Indeed, an examination of 
the assumptions carried out to arrive at case (c) one may note that A << 0.1 or 
at least of the order of 0.20 for the relation to hold reasonably well. Choosing 
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a value of A = 0.2, we note that, 

for X/A = l, (mS = 7.75, A = 0.2) R x = 1.4xl0 5 (2.19) 

for X/A = 2, (m6 = 9.8, A = 0.2) R x =8.3xl0 5 (2.20 

(A choice of A - 0.2 already introduces an error of about 8% in the approxi- 
mation of the boundary condition (4.8) in Benjamin's paper.) 

A choice of X/A - 1 implies a boundary layer growth of 40$ (X/A - 2 
implies a 22.5% change in 6 *) thus violating the assumption of parallel flow. 
According to our numerical experiments a 40% change in 6* introduces a change 
of similar order in the value of pressure due to nonparal lei ness of the flow. 
Thus for case (c) to be valid (X/A > 2) R^ has to be an unacceptably high value 
(of the order of millions). 
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Comparison of Navier Stokes Solutions with Benjamin Theory 

In order to verify the appropriateness of Benjamin theory, computer 
simulations were conducted using the Navier Stokes solver. The results of 
these calculations are summarized in Tables 1-1 through 1.3 (for a fixed x/A 
- 10 ). 

Table 1.1 presents results for a starting Reynolds number R x = 100; 
for this case 6 /A - 5.0. Results are presented for sinusoidal wave amplitudes of 

10- 4 , 10- 3 , lO" 2 , 2 x 10- 2 , 5 x 10- 2 , and 10" 1 . 

Comparison with Benjamin's theory indicates that both the linear profile 
and the sinusoidal profile approximations show poor agreement with the 

present computational solution for pressure and shear stress, even for 

small amplitude waves. 

In Table 1.2 we give the results of computer simulations at = 1000 
with X/ X = 10. Here again Benjamin's theory and the present Nav ier-Stokes 
results exhibit about a 30% difference in. pressure amplitude. The 
sinusoidal profile approximation appears to agree better than the 
assumptionof a I i near prof i I e in terms of pressure phase. 

In Table 1.3 results are given for R = I0 4 . Here the agreement 

X 

between theory and simulation is excellent for C^. The linear profile 
assumption seems to predict the shear stress distribution better than 
the sinusoidal profile approximation. This agreement may be fortuituous. 



A further aspect that needs to be checked regarding these comparisons is 
the role of (SA) . Table 1.4 Indicates comparisons for various 6/X's 
to obtain a reasonably large 5, we have chosen hypothetical values of 
R = 10 4 , U = !0 2 In/sec x = 100 in, etc In the simulations. We have 

4 

performed two sets of Navier-'Stokes calculations here. For 55 10 a 

Blasius profile Is used as the starting profile In the first set o f calculation 

and for the second set the oncoming flow Is taken as a shear flow with the 

same U’ (0) at the wall as .the Blasius profile. [In the freestream I.e., 

y> !/k (where k = UMOl/U^ ) the velocity Is kept ~ U , and f6r y<J/K,u =u ! (0)y.] 

It is seen that the Navier Stokes solutions differ from each other by 
10%. This is to be expected due to the difference of profile shapes. 

However Benjamin's theoretical predictions do not agree with 
the Navier Stokes solution. For 6 /?, = 1.0 both the sine and linear 
approximation agree within 10% in amplitude. For larger values of 6 /x the 
discrepancy is greater. 

We conclude from these numerical studies that the asymptotic 
theoretical results are evidently not accurate for moderate 
Re. 

Consideration of Couette Flow Assumption for Navier Stokes Solvers 
A number of works (theoretical methods and Navler-Stokes 
solvers) (Refs. 9,10, ID utilize a shear flow assumption 
for the unperturbed flow. These computer algorithms 

use periodic inflow and outflow condition with a moving lid at the top. For 
this case it is well known that a true steady state is available. 

It is, therefore, interesting to examine how well these codes may 
perform in the prediction of developing flow over a wavy wall. 
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Since the flow is assumed to be a high Reynolds number flow, there should 

be no dependence on development length (X q ). Thus the 
characteristic Reynolds number is R A = U a a /v where is chosen as a 
velocity scale at a distance a from the wall (U A = kA where k is the 
velocity gradient at the wall). We observe from Table 1.4 that a 
linear profile for the primary flow produces an error of about 10% In 
results when compared to the Blasius profile whenever the wall velocity 
gradient is chosen the same for both cases (see the Navier Stokes results of 
Table 1.4). Since the wall velocity gradient of a laminar flow at a Reynolds 
number R x = IL X/v , is given as 

U‘(0) = 0.332 It/Uo/vx = 0.332 lk/x Rl/2 (2-21) 

this must be equated to < in order to make an appropriate couette flow 
approximation of the developing boundary layer. 

Thus K = 0.332 IL/x Rl/2 (2.22) 

and R A = kA 2 /v = 0.332 lU 2 /vx R+ 1 / 2 (2.23) 

= 0.332 R3/2 (x/a J- 2 

Therefore a Couette flow Reynolds number corresponds to a Reynolds number 

R x = 2.0856 l\ 2 / 3 (x/A ) 4 / 3 (2.24) 

In Table 1.5, for a given R A = 5000 and waveheight a /A = 10“ 3 »we have 
performed calculation of Cp o and ^ f or various (x/a ) 1 s . [For R A = 5000 
and x/A = 1 Eqn. (2.24) gives R x = 609.8; for x/a = 2.0 R x = 1537; and 
for x/a = 4, R x = 3872.] 
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If the couette flow results are to be meaningful, the Cp values 
should agree with one another. In reality they vary widely depending on 
the value of R x or x/x . The Implication Is that It Is not the 

shear flow assumption itself that is at fault, but the denial of the memory 
to the fluid (development length) by virtue of defining . 
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Influence of Surface Wavi nes s on Flow Stability 


Laminar flows are subject to instabilities and eventual 

transition to turbulence at large Reynolds numbers. At very low 
Reynolds numbers disturbances introduced into the flow decay with 

time. When the Reynolds numbers are above the critical Reynolds number 
(R^ * = 520 on a flat plate) small amplitude disturbances may grow in time 
depending on the nature of the disturbance (wavenumber or frequency). This 
is an eigenvalue problem i.e. 

F(o^ , cu) = 0 (3.1) 

where tX is the wavenumber and w is the angular frequency (complex number). 
The initial growth of a small amplitude disturbance can be studied using the 
parallel flow assumption for the primary flow. The equation system governing 
the disturbances is the classical Orr-Sommerf eld equation 

(D 2 - a 2 ) 2 $ = iotR [(U-c) (d 2 - a 2 ) $ - u"$] (3.3) 

In general. Squire’s theorem states that three dimensional disturbances 
are typically less disruptive than the two dimensional disturbances 
small amplitude fluctuation in incompressible flow. 

The linear growth regime of small distrubances had been studied 
exhaustively (6) for flat plates. In general, the source of the 
initial disturbance is very difficult to control even under rigorous 
laboratory conditions. Abrupt or drastic changes in surface 
geometry (trips, roughness) sound sources in the field or even 
freestream vorticity can provide the initial disturbance which 
causes the flow to undergo transition. 


16 





We next examine the subsequent distrubance growth behavior 

for a disturbance free stable laminar flow on a flat plate 

encountering a wavy wall of a given amplitude and wavelength. 

In section (2.4) we stated that the solution of the Navier- 

Stokes equations was obtained for flow over a wavy surface with 

a starting velocity profile, i.e. a Blasius profile at time 

t = t , 
o 


i .e . , 


U(X,Y,t 0 ) = U 0 (Y,t 0 ) 


(3.4) 


Once the calculations begin, a perturbation field is generated 

due to the change in geometry. For instance, the perturbation 

to U at (t + At) is at most of the order (kaAt) where the 
o o '■ 

wall has a wavenumber f k’ = 2ir/X. For the next few time units 
the flow has to undergo an adjustment towards the eventual 
velocity profile U(X,Y,t) appropriate for the wavy surface. 

We note that Benjamin , using asymptotic theory, has shown that 
the perturbations from a mean flow over a wavy surface with 
a « 1 , and ka << 1 satisfies the Orr-Sommerf eld type 
equation ((see 3.1) of Ref. 3). The solution of the Navier- 
Stokes equation at early times can therefore be viewed as 
perturbations with the wavenumber *k’. Therefore it is reasonable 
to expect that the solution at initial times of the unsteady 
Navier-STokes solver can be studied to get information regarding 
the stability of the flow. 
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Results of Computer Experiments 


When the Reynolds number of the oncoming flow is well below a critical 
value the presence of waviness does not affect the stability of the flow. As 
an illustration the plot of pressure at a point with time (non-dimensional 
with x/2irlto) is shown in Fig. 3.1. The calculations were performed under 
the following conditions: X 0 = 1000, x = 12, d* = 10, a/x = 10“ d . The 

initial adjustment time from the figure is roughly 24 units. 

Since small amplitude perturbations of the Orr-Sommerf eld equations travel at 
a speed of roughly 0.2-0.25 U at these low Reynolds numbers, this time 
corresponds to a traversed distance of one wavelength for the perturbation. 
[24 x (. 25LL ) x (X /2n LU ) * Xj 

After a time of 24 units the pressure shows a slight decay with time. 
This decay has been found to be strictly porportional to the growth of the 
boundary layer in all of our calculations, i.e.,A 6*/6* * Ap/p. For the 
problem considered AS*/<$* = 1.8% and Ap/p - 4/438 = 1%. 

In Fig. 3.2 a plot of the pressure at a point is presented for the case 
R x = lCf 4 , X = 10^ , X = 60. Again one can observe that the initial transient 
settles down by a time t = 30 units. The decay of the pressure for 
subsequent times can again be compared to growth of the boundary layer. From 
the figure we deduce that the boundary layer grows by AS* = 8.6%S* and the 
pressure decays by Ap = 11. 7%p . 

When the Reynolds number of the oncoming flow is near the critical value 
the pressure transients show different patterns. Figs. 3.3 
through 3.9 show the behaviors of the pressure transient with time for an 
oncoming Reynolds number R x = id 5 (X 0 = 10" U* = 216"/sec.) and various 
wall wave combinations. 
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Figure 3.3 presents results for a wall valve of amplitude a = 5 x 10 


-4" 


with a wavelength A = 2" . The nondimensional wavenumber a= 


2tt<S* 


for this 


case is a =0.1709. The pressure transient contains a mean component which decays 
with time along with a Tollmein-Schlichting wave which behaves like a damped oscil- 
lator. Notice that the successive maxima are equally spaced in the time domain. 

The mean line dividing the envelope of maxima and minima for this case is the 
true mean pressure due to the wavy surface. The period of separation between 

successive maxima for this case is t_ = 18 units. [t =T A/2 tt U and hence 

0 o o °o 

2 

(jd = 2TT/t = 4 tt U /Ax .] Thus the nondimensional frequency 0) 6*/U = 0.0596 
r'o °°o ^ J r co 

give c r “ w /a = 0.349 U^, the velocity of propagation of the T-S wave. The 

logarithm of the ratio of successive maxima yields io.t . For this case 

- l o 

_3 

= -4. 4 x 10 . A separate Orr-Sommerf eld calculation for 

a linear perturbation to the flow at R^ = 10~* gives the value 


a) 6* = 0.0590 

r 

U„ 


U 3 6* = -4.226 x 10 
i 

U 


-3 


Figures 3.4 and 3.5 show plots of pressure transient at two points separated in 
space by a distance A/8, for the case a=l x 10 \ A=l" (R x = 10^, X = 10", = 216 M / 

sec). The nondimensional wavenumber a = 2 tt 5*/A = 0.34175 for this case. The suc- 
cessive maxima and minima are separated by a time period of 16 units. 

-3 

yielding (0^=533 rad/sec or 0)^6* /U^ = 0.1363; for this case (*) 6*/!!^ = -7 . 5 x 10 
The velocity of propagation 0 ^= 0.399 U m . By comparing Figs . 3.4 and 3.5 it is 
seen that the pressure signal in (3 . 5) is shifted by a time unit = 2 units from that 

in (3.4). Thus 

A/8 = c t giving c = tt/8U = 0.393U 

r j. r w oo 
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Solution of the Orr-Sonmierf eld equation yields 

^ 5 * /U oo = 0 * 13743 ; w.<$*/U =-2.1865 x 10“ 4 

r w X 00 

In fig. 3.6 a plot of pressure transient for a wall wave of amplitude 
a — x ^ 1 is given (R^ = 10 , X = 10", = 216"/sec) . From the figure 

uj = 523 rad/sec w 6*/U = 0.132 and w. 6*/U = -7 . 92 x 10" 3 and c = 0.386 U 

■*- r 1 00 00 

-2 

For Fig. 3.7 the amplitude of the wave wall was increased to a=l x 10 in. 

The separation between successive maxima no longer display constancy in period . 
even though they are approximately equal. Indeed the nonlinear nature of the 
wave is barely observable (shorter fetch during rise and longer fetch during fall 
in the oscillation in the initial phases of the oscillations) . 

In Fig. 3.8, the pressure transient for a wall wave of amplitude 
a =3 x 10 ^ in, A = 1" (R^ = 10 3 , X = 10", U = 216"/sec) is plotted. Notice the 
growth of the pressure perturbation with time. 

-2 

In Fig. 3.9 the pressure transient for a wall wave of a = 5 x 10 in 
(A = 1, R = 10 3 )indicates the essentially nonlinear nature of the growth of the dis- 
turbances. For both (3.8) and (3.9) stable quasilinear solutions are not avail- 
able indicating the time unsteady separation of the flow in the wall region [the 
solutions indicate the growing nature of the separation bubble eventually leading 
to complete breakdown of the flow at a later time] . 

In Figs. (3.10) through (3.13) the flow conditions are r =5 x 10 5 : U = 1080"/ 

x 00 

sec, X = 10". The wavelength of the wall wave is chosen as A = 1.0". 
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Fig (3.10) shows plots of pressure for a wave amplitude a=l x 10 . From 

the figure, a = 0.1528; a) <5*/U =0.04796; a).6*/U =-2.9 x 10 ^ and c =0.314. 

° ' »j*O0 'loo j- 

—4 

The Orr- Sommer f eld solutions yield U) <5*/U =0.046275; w.6*/U =-2.097 x 10 for 

j r 00 i oo 

this case. 

_3 

In Fig. (3.11) the amplitude of the wall wave is a =3 x 10 . At this 

amplitude itself the nonlinear effects are barely visible. 

-2 

Fig. (3.12) depicts the situation for the case of wave amplitude a = 1 x 10 
The wave form now consists of a group of traveling waves. This is ob- 
tainable by considering Fig. 3.12 along with Fig. 3. 13 where the pressure 
signal at a point A/8 downstream is plotted. The constant phase shift can be 
converted to obtain the group velocity C - 0.26 U , At higher amplitudes 
a total breakdown of the flow is noticable in convective time scale of 
0 ( 1 ). 
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Methodology for Evaluating the Drag Reduction Characteristics of Wavy 


Surfaces 

In a wind tunnel the drag reduction characteristics of various nonplanar 
surfaces may be compared by measuring the total drag forces on the surface 

and comparing these to a planar surface. When numerical simulations 
are used to evaluate the drag characteristic of various nonplanar surfaces, 
the following requirements must be satisfied: (1) The simulation should be 

highly accurate and capable of predicting the skin friction and pressure drag 
on the surface accurately. (2) The simulation should be more economical 
than the experimental work. 

The Navier-Stokes algorithm reported in (Refs. 1 and 2 ). 
satisfies both these criterions. The Fourier Spectral method has been found 
to predict the total drag on wavy surfaces in turbulent flow in agreement 
with experiments (Ref. 2). The algorithm reported in (Refs. 1 and 2) uses 
periodic boundary conditions in the flow direction. While this has a decided 
advantage in terms of cost effectiveness as compared to an inflow/outflow 
steady state solver (spatial), there is a nagging doubt with 
regard to the accuracy of the solution vis a vis the inflow-outflow 
codes. The results that are reported herein should provide some answers 
to these doubts. 

We note that the inflow-outflow codes can be a burdensome asset when flow 
over more than one wave (multiple waves) needs to be studied. For each wave 
a resolution of at least 30 points are needed in the wave-direction (in order 
to estimate phase of pressure distribution with an accuracy better than 5 U ) 
thereby increasing computation time as well as storage requirements. 
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Secondly, conformal mapping techniques become prohibitively expensive to 
employ for these problems because a flow over N waves increases the 
the computational labor N fold due to mapping complexity. A spatial code in 
these situations is less desirable than a properly constructed temporal 
code. The Navier Stokes solver described in Sec. 2 is a temporal 
code. In order to correlate the temporal evolution of solution to spatial 
evolution of a developing flow, we consider the following three 
illustrative examples: 

1 ) Consider at time t = t^, a flat plate which is set suddenly in motion 
with a velocity CJ = -U» . By Galilean transformation this is equivalent to 
the flow surging at U = at t = t Q and the plate stationary at 
subsequent times (a classical Rayleigh problem.) For this problem the 
governing equation of motion is 

9 u/9 t = v du Z /2y Z (4.1) 

with, at t = 0: U = Uo v > 0 

U = 0 at y = 0 

and U = 0 at y = 0 at all times. 

The solution for velocity, shear stress, etc. is 

U = llo erf y// 4v"t (4.2) 

and T w = Uo/pp A t (4.3) 

where erf(z) = 2 // it e“^ (4.4) 
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A Reynolds number can be constructed for the flow as, 

Rr = U* 2 t/4v (4.5) 

The displacement thickness d x for the flow is 

6i - 1.99 (Uot)//^ (4.6) 

and the skin friction 

C f = 0.5624//R7 (4.7) 

For a developing Bias i us flow the relations are 

= 1.72 (4.8) 

and Cf = 0.664//R x (4.9) 

Because of the difference in the numerical constants in (4.5) throuqh 

(4.9) a proportionality between X and t is not strictly correct. (An error 

of 10% occurs for skin friction values if the relation X = Ubt/2 = t 

is used for Rayleigh analogy. ) If, hcwever, a velocity scale is chosen such 

that Cq - .664/. 5624 Cq - .6U» the Rayleigh solutions give reasonably 

accurate values for various profile parameters. Thus even for this problem a 

meaningful analogy can be constructed. 

In Fig. 4.1, a plot of [(C f /C f ) z - 1] vs T is shown. The solutions 

o 

were obtained using the periodic Navier Stokes Solver with amplitude of wave 

a = 0 and with the intial condition given by (4.1).* Our calculations qive 
Cq - 0.620, , in good agreement with the exact solution. 


*The group velocity of computation can be determined by using the 
formula, = [0.4409 Kv/ ( where K is the slope of the 

curve shown in Fig. 4.1. 
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To investigate the nature of the solution with time the following 
initial conditions are prescribed 

t = tQi u = U(Y) , where U(Y) is the Blasius profile corresponding to sane 
x Reynolds number Rj^* 

Since the equation system (4.1) does not have an x dependence at t = t Q , 
the solution at all times will also be x- independent . Because of the parabolicity 
of (4.1) at large times the solution tends to the Rayleigh solution. Of 

interest, however, is the nature of the solution at small and intermediate 
times. In the small to intermediate times memory of the initial 
profile exists for the flow. 

We illustrate this behavior for the solution through Fiq. (4.2) and 
(4.3 ) s where R x 0 °f 200, 1000 are chosen in the examples. The group velocity 
of computation C is equal to 0.21, which is the 

(j 

perturbation speed of propagation at these Reynolds numbers.'* The time domain 
in Figs. (4.2) and (4.3) is normally the time taken to obtain a quasisteady 
solution for a wavy wall of the same periodic length as the computational 
box. 

There is hardly any scatter of data points from the straight line in 
these figures. This means that the solution of the flat plate at different 
times (t Q + t x ) can be related to the solution for a developing boundary 
layer at Xq + x, where x L = . This concept is of crucial interest 

in developing a raethology for rating the drag characteristic of various 
nonplanar surfaces. 

*The group velocity of compu t at jpn can be determined by using the 
formula, C Q = [0.88l8 Kv where K is the slope of the 
curves showrf in Figs. 4. 1-4.3. 
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The approach that we have adopted for draq ocarparisons is to start a wavy 
wall calculation at time T 0 with a Blasius profile corresponding to an X 
Reynolds number RXq- At various time intervals (t 0 + r&t) the total drag 

of the wavy surface is computed. These quasisteady values are compared to 
the local skin friction coefficient value at (X Q + nCgAt) in order to 
obtain a reasonable comparison of the performance of various surfaces. 
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Nonplanar Surfaces for Drag Reduction in Laminar Boundary Layers 


Turbulent flow covers much of the surface area on CTOL aircraft. 

However, the residual areas of laminar flow (e.g., near nacelle leading edges 

and the fuselage nose region) can contribute large skin friction drag levels 

due to their initially low Reynolds number. Research reported herein 

indicates that certain nonplanar geometries can reduce net skin 

friction drag, even in laminar flows. Thus nonplanar waveforms may 

be a viable alternative to the more esoteric technique of 

laminar drag reduction through "supersmooth" ("slip" boundary) 

surfaces. A computer simulation technique (Navier-Stokes Solver) 

using high order numerical methods has already been developed for 

both laminar and turbulent flow wavy wall cases (Ref. 1,2) and 

it was found that the computer simulations accurately predict 

(for turbulent flows), within a 4 percent error margin, drag 

levels measured in the wind tunnel. 


In the laminar flow region the prediction capability of the present 
code is further improved due to the obvious absence of turbulence modeling and 
associated inaccuracies. Fig. 5.1 indicates four waveforms with a/X = 0.2. 
Calculations of total drag for these waveforms were conducted in the Reynolds 
number range R x - 100 to 1000. The variation of total drag as compared to 
flat plate drag is shown in Fig. 5.2. The unsymmetric V-groove of Figure 
5.1a was found to have a net 17 percent drag reduction compared to a smooth 
flat surface with the same planform area. The sine waveform was shown to 
have a drag penalty starting with the lowest Reynolds number. 
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Our conclusions are that, at least at low Reynolds number, certain 
nonplanar asymmetric waveforms may perform better than a smooth flat plate. 
The mechanism of viscous drag reduction is the nonlinear effect of the 
pressure distribution over these waveforms. The asymmetric shape of the 
surface is curcial in providing a net drag reduciton. Ony in certain cases 
(Figs. 5.1a,b ) is the associated pressure drag on the wave smaller than 

the skin friction reductions. These studies are ongoing, with emphasis upon 
higher Reynolds number and further surface optimization. 


28 


The Effect of Suction and Injection on the Total Drag of a Wavy Surfaces 


The increased cost of fuel has reemphasized viscous drag reduction 
research. The skin friction reduction obtained through surface mass 
injection is well known. However, the effect of spatially non-uniform 
surface injection on the amount of drag reduction (with fixed mass flow rate) 
has not been examined in detail. An interesting approach to providing local 
(passive) mass sources and sinks is the use of a porous wavy-wall 
configuration with its alternating pressure levels to provide self -induced 
mass transfer through the porous wall. In this case the pressure 
distribution over the surface wave can result in a non-negligible pressure 
drag which may be a significant part of the total drag on the surface. The 
purpose of the present study was to examine theoretically the effects of 
various spatial variations of suction and injection on the drag of both flat 
and wavy surfaces. 

Calculations have been made for two (four foot long) surfaces, one flat 

and the other sinusoidal with a wave length of one foot. In each case the 

undisturbed velocity was 50 fps, the initial boundary layer was turbulent 

at Rq = 4500, and the initial boundary layer thickness was approximately 

1.8 inches. These conditions correspond to the experiment of Sigal for a 

14 

wavy-wall with no suction or injection . For the flat surface, a 
finite-difference boundary layer code described in Ref 0-5) was used to 
determine if spatial variations in mass transfer for a fixed total flow rate 
could result in enhanced drag reduction. This code was used previously^ to 
study the effect of Mach number on drag reduction due to uniform surface 
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injection. For the sinusoidal surface, calculations were made using the 
Navier-Stokes code for a sinusoidal distribution of suction and injection 
perpendicular to the wall. In both programs modelling constants of k = .41 
and (£/6) max = .085 were used in zero order (mixing length) turbulence 
closure schemes, and the wall damping factor, A + , was assumed to be the 
following function of V -1- for low speed flow : 

A + = 26 exp (-5. 9 V + ) (1) 

8 13 . . 

A number of works » have suggested using A" r = fCP" 1 ") and k = f(P) 

corrections as a means of including the equilibrium effects of pressure 

gradient on turbulence. However, computations reported in Ref. 2 indicated 

that the pressure and skin friction distributions over Sigal's wave could be 

adequately predicted using zero pressure gradient constants (probably due to 

the nearly "frozen" nature of the turbulence structure); therefore, no 

corrections of the type proposed in Refs. 8 and 13 were applied in the present 

calculations. 

The Navier-Stokes code was tested by comparing calculations on a flat 
surface with a sinusoidal distribution of injection and suction (v w /uoo = 

.005 at maximum suction and injection) to similar calculations made with the 
boundary layer code. Results obtained for an impermeable wall, a flow with 
constant blowing, and the sinusoidal suction solution for one cycle of 
suction and blowing are shown in Fig. 6.1. The drag coefficient calculated 
for the surface element (four cycles of suction and injection) was 2.86 x 
10~~* from the boundary layer code and 2.81 x 10~~* from the Navier-Stokes 
code. Note that the integrated skin fricion drag over the region where the 
sinusoidal suction was applied is greater than the drag of an impermeable 
surface. 
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In Fig. 6.2, the Navier-Stokes calculation for a sinusoidal wave (h = 
0.01“, A = 12", 6 =3") with suction velocity (v w /u 00 ) max = .003 is 
presented. From the figure it is apparent that the self bleed situation (6 = 
0°) imposed a substantial drag penalty over that of a wavy surface with no 
mass transfer, mainly through an increase in pressure drag (phase shift). 

The opposite situation occurs for peak suction at the crest (<Ji = 180°). Here 

/ 

the reduction is mainly due to a pressure thrust. Calculations were also 
made with different peak suction rates for the same wave. All these 
calculations indicate similar trends with $. However, there seems to be a 
critical (v w /u„ ) max for a given wave and flow conditions for which the 
total drag reduction is a maximum when ® 180°. Above this critical value 
the total drag goes up for all $ • It may thus be possible to obtain a total 
drag reduction over the wavy surface by a suction distribution with 6 - 
180°. However, the associated penalty for collecting and ducting the flow 
(which is not considered in this analysis )implies that this is not a 
desirable option for drag reduction on sine waves. The original possibility 
which prompted the study, that of using passive bleed over the wave, is seen 
to be detrimental. 
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Drag Characteristics of Nonplanar Waveforms in Turbul ent Flow 

The Navier Stokes solver discussed in Ref. 1 and 2 has been used to 

study the drag characteristics of various nonplanar waveforms in turbulent 
flow. For turbulent flows closure of the time average Navier Stokes 

equations is required. In Ref. 2, various closure models have been tested. 

It may be recalled here, that, based on the available experimental data, the 
zero equation model with equilibrium constants performs adequately in flow 
predictions for nondimensional amplitudes of ka up to 0.2. [The experiments 
against which companions were made in Ref. 2 are (i) Kendall’s experiment (ka 
= 0.19, Rq - 2000), (i i ) Si gal's experiment (ka = 0.17, Rq = 6000) (i i i ) 

LaRC inhouse experiments in the 7 x 11 tunnel (f^ = 2000, ka = 0.03- 

0. 125, with sine waveforms and other nonplanar waveforms. )J The design 
studies that are given here, therefore, use the zero equation closure model 
((£/<S) m ax = 0.09; K = 0.4; A + = 26.0). The methodology discussed in Sec. 

4, is used to determine the drag characteri Stic of the nonplanar waveforms: 

1. e., the evolution of the solution over the nonplanar waveform is compared 
to the evolution over a flat surface. The design studies are purposely made 
for the 7" x 11" tunnel conditions of LaRC so that "promising" candidate 
surfaces can be tested at a later date. (These experiments are currently 
underway . ) 

Based on the results of earlier experiments it was observed that 
symmetric sinusoidal waves have higher form drag penalties compared to other 

waveforms of the same wavelength and amplitude. The purpose of this study 
was to compute various asymmetrical surface shapes and study their overall 

drag characteristics. j n following we report the results of our 

investigation. It was noted in Sec. 4 that the convection speed of the 
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calculations (i.e., Cq) for turbulent flows was approximately Cq ~ 

0.8LL. The initial (startup) transients decay at a time scale 
smaller than for laminar flows. In fact, we note that the skin fricion and 
pressure distributions settle down after 4 units of non-dimensional time 
(4(x/2TTUx>)sec) i.e., in a time required for the box to float a distance ~ 

0. 5x . 

In the following we briefly describe the surface waveforms drag 
characteristics. 

(1) Initial (Hat Plate) Results 

The results of computer solutions for a flat plate at 7 x 11 conditions 
(IL = 75"/sec, <5 initial = 0.375") at various time evolutions 
(non-dimensional time) are presented below. 
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Entry B1 


Surface: sine wave, a = 0.03", X = 1" (ka = 0.188) 

Results of the simulations are 


Time 

3.6 

4.2 

4.8 

5.4 

6.0 

C n x 10? 
u f . 

3.645 

3.585 

3.528 

3.475 

3.428 

Cn X 10* 
P 

2.209 

2.027 

1.985 

1.988 

1.996 

C Dj x 1C? 

5.851 

5.613 

5.513 

5.463 

5.424 

% 

Drag Change 

59.2% 

54. 1% 

52. 5% 

52.3% 

52.3% 


There is a drag increase of 52% for this wave which is mainly due to the 
pressure drag on the surface (form drag - 56% of the flat plate). 
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Entry B2 

Surface: sine wave, a = 0.015", X = 1" (ka =: 0.10) 

Results are tabulated below 



% 15.35 14.0 13.4 13.2 13.0 

Drag Change 
































Entry B 4 

Surface: symmetric V groove, a = 0.0236, A = 1" 

(The V-groove has the same wall slope as the maximum slope on a sine surface 
(i.e., Case B2) for this case.) 


Drag Change 


3.6 


i 

4.8 

5.4 

3.600 

3.53 

5 

3.472 

3.413 

1.077 

1.00 

5 

0.979 

0.974 

4.677 

4.54 


4.45 

4.387 

27.3 

24.6 


23.2 

22.3 


Note that the pressure drag is proportional to a 2 for a given k, 
simulation (compare B3 and B4; B1 and B2) attest to this. 
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Entry B5 

Surface: unsymmetric V-groove, a = 0.015", X = 1" 

The unsymmetric V groove can be analytically described as; 

y = a(l - 2.857 x/x ) 0 < x/x <. 0.7 
y = a(6.67 xA -5.67) 0.7 < x/x _< 1.0 


Time 

3.6 


4.8 

5.4 

6.0 


3.649 

3.604 

3.561 

3.521 

3.484 

HBBI 

0.446 

0.409 

0.388 

0.376 

0.366 

team 

4.096 

4.013 

. 

3.950 

3.897 

3.850 

% 

Drag Change 

11.45 

10.20 

9.30 

8.64 

8.05 



















Entry B6 

Surface: unsymmetric V-groove, a = 0.03", x = 2" 

The surface has the same shape as Entry B5 


Ti me 

3.6 

4.2 

4.8 

5.4 

6.0 

C n x lti* 

_f 

3.499 

3.455 

3.418 

3.386 

3.359 











Entry B7 

Surface: symmetric rounded and straight shape. 

Ihe surface can be described best as follows. The crests have a convex 
curvature k = 0.5/in. and the troughs have concave curvature k = 1.0/in 
rounded protions are connected by straight line segments; a = 0.015" * : 


Drag Change 


5. 


3.552 


0.470 


4.022 


12.13 


m 


6 . 


9 


68 


3.986 


1 
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Entry B8 


Surface shape same as B7 except that convex curvature is k = 1/in. and the 
concave curvature « = 0.5/in. (a = 0.015" A = 1") 











Entry B9 


Unsymmetric Rounded shape. 

This surface has a convex curvature (at the crests) k = 1/in. and concave 
curvature k = 0.5/in. Crest to trough distance in downwind direction is 
0.75 a. The curvature extends to O.lx on either side of crest and trough 
followed fay straight line segment (k = 0). The extent of straight line 
segment is 0.6a from crest to trough and 0.2a from trough to crest. 


Time 

3.6 

4.2 

4.8 

5.4 

6.0 

taum 

3.474 

3.408 

3.326 

3.250 

3.180 

IMfl 

0.453 

0.383 

0.344 

0.324 

0.314 

HB 

3.946 

3.790 

3.670 

3.574 

3.494 

% 

Drag Change 

7.37 

4.03 

1.55 

-0.4 

-2 


Ihere is a possibility of drag reduction for this case! 





















Entry BIO 



The surface shape is similar to B9 except that the curvature on the crest 
side is now 0.5/in. and at the trough side is 1/in. 
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Entry Bll through B13 are sinusoidal surfaces that are asymmetric with a 
0.015", x = 1" 


Entry Bll 
Surface 


y - a cos (k^x) 0 < x/A < 0.75 k A = 2/ 3k 
y = a cos (k z x) 0.75 < x/A _< 1.0 k z - 2k. 


Time 

3.6 

4.2 

4.8 

5.4 

6.0 

C n x 10* 
J f 

3.629 

3.580 

3.534 

3.492 

3.453 

C D x 10* 

0.522 

0.451 

0.412 

0.388 

0.372 

C n ' x ltf 
U T 

4.151 

4.031 

3.945 

3.88 

3.825 

% 

Drag Change 

13.0 

10.7 

9.16 

8.17 

7.35 



Entry B12 


The surface shape is the same with the flow direction reversed. 
1 i.e. y - a cos k z x 0 < x/x _< .75 

y = a cos k i (x - . 25X) .75 < x/x _< 1.0 

F 


Time 

3.6 

sa 

4.8 

5.4 

6.0 

C n x ltf 
u f 

3.675 

3.632 

3.595 

3.563 

3.533 

C D x ltf 

-P- 

0.993 

0.979 

0.998 

1.010 

1.023 

V x lo- 

4.668 

4.611 

4.593 

4.573 

4.557 

% 

Drag Change 

27.1 

26.6 

27.1 

27.5 

27.8 


Notice the drastic increase in total drag compared to Bll. 
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Entry B13 


I he surface can be described as 

y = a cos kx 0 _< x/A <_ 0. 5 
y = -a(l+2coskx/2) 0.5 _< x/X <_ 1.0 


Time 

3.6 


4.8 

5.4 

6.0 

C n x 1C? 
u f 

3.642 

3.598 

3.557 

3.518 

3.482 

Cn X 1C? 
P 

.542 

0.503 

0.490 

0.487 

0.487 

Cn' X 1C? 
U T 

4.184 

4.101 

4.047 

4.005 

3.969 

% 

Drag Change 

13.9 

12.6 

12.0 

11.7 

11.4 




















Entry B14 through B17 are hybrid sine straight line surfaces 


Entry B14 


a = 0.015" x = 1" 

y- a(l - 2.66x/x) 0 < x/x < 0.75 

y = a cos (4 tt(x-0. 5) /X) 0.75 < x7>. < 1.0 






















Entry 1315 


a = 0.015" x = 1 


Surface: y - a (a - 4x/x ) 0 < x/x <0.5 

y = a cos 3rx/x 0.5 < x/x <1.0 


Time 

3.6 


4.8 

wwm 

6.0 

EHEB 

■Ufa 

3.667 

3.628 

3.591 


3.524 

UNKifl 

0.484 

0.442 

0.424 

i 

0.408 

■HUH 

4.151 

4.070 

4.015 


3.932 

% 

12.99 

11.7 

11.1 




% 

Drag Change 


12.99 


11.7 


11.1 













































Entry B17 

y = 
y = 

a(l - 13.3X/A ) 0 < xA < 0.15 

a cos (x + 0.7Q\ ) k x 0.15 < xA < 1.0 

where 

ki - 

( 1 . 176rr A ) 


Time 

3.6 

KOI 

4.8 

5. 

4 

6.0 

ESBH 

3.878 

3.898 

3.917 

3.93 

3.940 

|H| 

1.048 

1.036 

1.051 

1.08 

1.145 

isom 

4.927 

4.934 

4.967 

5.011 

__ 

5.054 

% 

Drag Change 

34.1 

35.4 

37.4 

39. 

7 

41.8 























Entry B18 



i 


* 

i 


Surface: y = a cos 2 ttx/a 

y - a(4x/x - 3) 


0 < x/x <0.5 
0.5 < x/x < 1.0 


Time 

3.6 

4.2 

4.8 

5.4 

6.0 

Cn X 10? 
u f . 

3.654 

3.612 

3.572 

3.535 

3.499 

C D x 10* 

fL .. 

0.542 

0.508 

0.496 

0.495 

0.496 

C n ' x 1C? 
U T 

4.195 

4.120 

4.096 

4.030 

3.995 

% 

Drag Change 

14.18 

13.1 

12.6 

12.4 

12.1 


Notice the slightly larger pressure drag for this case 
the flow encounters the straight segment first and the 
1 ater. 


compared to B15 where 
sinusoidal region 
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Secondary Instability of Wall Bounded Shear Flows 


The process by which a laminar shear flow undergoes transition to 
turbulence is of major interest in fluid dynamics. Here we study the initial 
transitional process that leads to eventual breakdown of wall bounded shear 
flows . 

The process by which a laminar flow may become turbulent, i.e., the 
transition process, depends greatly on the geometry, initial conditions and 
other specific features of the flow. For flows like Benard Convection and 

Taylor couette flow a sequence of ordered states of transition may be 
discerned. With increasing Reynolds numbers, chaos results in these flows 
through either (i) infinitely many period-doubling bifurcations; (ii) inter- 
mittent chaotic states; or (iii) three incommensurate Hopf 

bifurcations. The choice of which route the flow undertakes in the 
transition process seems to depend on the specific details of the flow, the 
geometry and initial conditions. This sensitivity is probably due to the 
comparati vely low Re numbers at which weak chaos appear in these flows. 

The situation in wall bounded shear flows is different. The transition 
to turbulence is violent in the sense that the flow snaps from a laminar 
state to a strongly chaotic state with no apparent intervening stable 
bifurcations from the mean flow. 

The instability proces can be categorized into 

(i) Primary (linear instability) of the basic shear flow 

(ii) Nonlinear saturation of the primary instabiltiy and the eventual 
formation of the secondary flow. 

and (iii) the secondary instability i.e. linear instability of the secondary 
flow. 
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A rational analysis of the physics of the instability requires 


understanding the roles of each of these three stages in the growth of 
disturbances that leads to breakdown of flow. In this section a summary of 
the work in progress in the area of stage ( i i i ) is presented. 

For analysis the primary instability was chosen as two-dimensional. The 
two dimensional primary instability results in a periodic nonlinear flow 
whereas the secondary instability likely leads to chaos. The secondary 
instability is assumed three-dimensional. The choice of two dimensionality 
of the primary instability is assumed with the understanding that though 
calculations of three-dimensional primary instabilities show an eventual 
transition to a chaotic state, the time and space scale of these 
instabilities do not bear connection to experimental 

observations. On the other hand the secondary instability considered here are 
inherently three dimensional, and have explosive growth rates (convective) as 
against primary instability (which grows in diffusive time scales). 

Method 

To study nonlinear saturation of the primary two dimensional instability 
in plane poiseuille flow the velocity is written of the form 

- N + 

v( 2 ) = \\ (z) x + 7 V n (2) (z,T),e ia n x 

-N 

The two dimensional velocity field may be specified as 

- N + 

1/(2) = (z) x + / v(2) (z,t ) e^ a ( x " c ^) 

-N 
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where the time scale t is for the propagation of the nonlinear waves and the 
time scale t is related to the attenuation of the travelling waves i.e. 

t = t/R (t is the diffusive time scale). The presence of the slow time scale 
in the evolution of general nonlinear two dimensional disturbances can be 
studied by perturbing the equilibria that are independent of x. In plane 

Poiseuille flow there are stable quasi-equilibrium states at sufficiently low 

Reynolds numbers i.e., R < 2900 and equilibrium states between 2900 < R < 

5000 suggesting the role of the slow time scale in this problem. On the 

other hand axisymmetric pipe flow is stable to all finite amplitude 

perturbat i ons . Our calculations suggest that for pipe flows the 
diffusive scale T does not play a role in the evolution of axisymmetric 
perturbation of the pipe flow. 

Results for Plane Poiseulle Flow 

The flows to be studied are assumed to be of the form of a general two 
dimensional velocity for which an infinite small three dimensional 
disturbance is superposed. 

i .e. , 

+ N + 

V = U (f (z) x + }' V n (2) (z,t) e la n x + cV( 3 ), e «1 

-N 

> +N . 

v(3) = )' V V nm (3) (z,t) e^n>y eloc n x 

m = ±1 n = -N 

i.e., the choice is of only one mode in the spanwise direction (due to linear 
separability). N modes are kept in the stream direction due to nonlinear 
effects. 
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Using the multiple (2 time) scales we write 


V = U|| (z) x + )’ V n ( 2 ) (z,x) e ia n( x-ct ) + 

n = -N 

N N 

e(e at v V Vnm( 3 ) (z,x) e 1 ^ e i «n( x ‘ ct )> 

m = ± 1 n = -N 

Ihus a is the growth rate of then three dimensional secondary 
instability. 

The results of our calculations are briefly summarized a follows for 
plane Poiseuille flow. For Re = 2900, for 3 greater than a threshold value ( 3 

« a for maximum growth) the three dimensional perturbation grows in a time 
scale (a - 0(1)). 

For analysis of the secondary instability mechanisms at lower Reynolds 
numbers > a primary mode at a high Reynolds number is chosen (e.g. R = 4000). 

It is seen that the three dimensional disturbances turn off below R < 400: as 
R -*■ oo the growth rate is independant of R and hence the instability is 
inviscid in nature. In particular for R* 1000 or less, viscosity significantly 
limits the growth rate. 
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Secondary Instabilities of Free Shear Flows 


The mechanism of secondary instability of free shear flows was 
investigated using the techniques discussed in the last section. The results 

for incompressible flows indicate that the secondary instabilities are 
convective in character, with a convective growth rate of 0.2. This is in 
agreement with Nishioka's experiments. The effect of compressibility on the 
secondary instabilities is also under current investigation. Work is also 
under progress on the Cyber 203 on compressible turbulent boundary layers. 
Black Box Version of Compressible Sally Code 

In recent years there has been attempts to design laminar flow control 
vehicles as a means of reducing fuel consumption of CTOL aircraft. The LFC 
technology has prompted renewed interest in the study of transition of both 
compressible and incompressible laminar boundary layers. The SALLY code was 
expressly developed as a design tool to perform efficient design calculations 
(eigen value calculations) of the stability of laminar boundary layers. The 
original version of the Sally code analyzes-, incompressible flows where the 
eigenvalue problem is the fourth order Orr -Sommerfeld equation system. For 
compressible boundary layers the eigen value problem is eighth order for a 
general three dimensional boundary layer. When the flow is two dimensional 
or in the absence of dissipation the eight order system reduces to sixth 
order. 

The basic equations for the linear stability analysis of parallel flow 
compressible boundary layers are derived using the small disturbance theory. 
A set of five O.D.E.'s are obtained under the assumption of locally parallel 
flow (i.e. 3 second order momentum equations, 1 second order energy equation, 
and one first order continuity equations). This system of equations can be 
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reduced to a set of eight first order equations making the system aminableto 
an intial value numerical integration. 

However, previous work in this area (i.e., initial value approach) 
suffers mainly in that often they are computationally slow and they require a 
reasonably good initial guess. Thus, they are not attractive as a design 
tool . 

In the present work the original system of five equations are solved 
using finite difference matrix methods. The norm momentum equation is not 
converted to a first order equation for pressure since this generates 
unstable spurious modes when the problem is done using a global method. 

Two eigenvalue search procedures are provided - (i) Global, when no 
initial guesses are available (i i ) local, when a reasonable guess is 
available and refinements are required. 

In the global approach care has to be taken to avoid the generation of 
spurious growing unstable modes that are physically not relevant. The 

spurious modes can be distinguished by refining the spatial discretisation; 
the number of true modes increase while the spurious modes remain the same. 

A key concept that was used to eliminate the generation of spurious 
results in this work is that one can eliminate the spurious mode by 
using a numerical method for the stability problem that gives numerical 
stability for the initial value problem also. In this situation no 
spurious modes appear. Details of this work are reported in Ref. (17). 
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TARLE 1.1 Rv = 100 
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TIME 











Fig. 3.4- Plot of pressure at point 1 on a wavy surface. 

a/A=10~ 3 ; x/A = 10; a = 0.34175. 
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TIME 







Fig. 3.7- Plot of pressure vs time at a point for a wavy surface 
R =10 5 ; aM = 10 _2 ; a = 0.34175; xM = 10. 






TIME 





TIME 









Fig. 3.11- Plot of pressure at a point on the wave, ; 

a/X = 3xlO~ 3 ; X/A = 10; a 0.1528. 
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Fig. 3.13- Plot of pressure at point 2 (A/8 downstream of point 1) 
R r = 5 x 10 5 ; a/A = 2 x 10“ 2 ; X/A = 10; a = 0.1528. 
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5.1- Plot of nonplanar geometries. 
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Flo S A Fig. 5. Id- Sine waveform 

18 * 5.1c- A symmetric v-groove. S 



Fig. 5.2- Plot of Drag characteristic of nonplanar waveforms in 
laminar flow as a function of the Reynolds number of 
oncoming flow. 



Fig. 6.1- Calibration plot of skin friction distribution over a wave cycle of suction and 
iniection distribution for a flow over a flat plate. 
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